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S: 

! In soft matter structure couples to flow and vice versa. Complementary to structural in- 

a ■ 

Q ' vestigations, we here are interested in the determination of particle velocities of charged 

o : 

colloidal suspensions of different structure under flow. In a combined effort of theory 
^ . and experiment we determine the Fourier transform of the super-heterodyne field auto- 

! correlation function (power spectrum) which in frequency space is found to be well sep- 

m ■ 

CN ' arated from homodyne contributions and low frequency noise. Under certain conditions 

^ . the power spectrum is dominated by incoherently scattered light, originating from the un- 

oo : 

Q , avoidable size polydispersity of colloidal particles. A simple approximate form for the low- 

' wavenumber self-intermediate scattering function is proposed, reminiscent to the case of 

j_j ■ non-interacting particles. We experimentally scrutinize the range of applicability of these 

simplified calculations on employing a parabolic electro-osmotic flow profile. Both for non- 
interacting and strongly interacting fluid particle systems, the spectra are well described as 
diffusion-broadened velocity distributions comprising an osmotic flow-averaged superpo- 
sition of Lorentzians at distinct locations. We discuss the performance and scope of this 
approach with particular focus on moderately strong interactions and on multiphase flow. In 
addition, we point to some remaining theoretical challenges in connection to the observed 
linear increase of the effective diffusion constant and the integrated spectral power with 
increasing electric field strength. 
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PACS numbers: 82.70.Dd,64.70.pv,81.10.Aj 



I. INTRODUCTION 



Soft condensed matter is typically characterized both by its softness and an internal structure 
on a mesoscopic scale. This allows convenient optical access to system structure and dynam- 
ics in equilibrium. Structural properties, obtained from static light scattering with good statistical 
accuracy can be compared to theoretical expectations and/or computer simulations |uil2Ll3|]. In par- 
ticular for systems of strongly interacting particles, local and global changes of structure and even 
phase transitions are frequently observed with and without flow, thus providing an interesting field 

nnnnnnn 

for studies of non-equilibrium phenomena [14J, |5|, la EL |8|, |9|, llOfl . We are intererested here with the 
complementary task of determining the (local) particle velocities in systems of interacting particles 
in arbitrary flow profiles (plate-plate or cone plate shear, Hagen-Poiseuille flow, electrophoretic- 
electro-osmotic flow in closed cells, thermophoresis, sedimentation etc.) and different structure 
(unstructured, fluid, crystalline) The combination of structure and flow measurements in suitable 
model systems should in the long run help to elucidate the fundamental coupling between both. 

A well developed arsenal of methods is available to study the flow behavior of fluids in general. 
Quite frequently tracer particles are used [jlll . Direct particle imaging nowadays is capable of 
fast scans of 3D velocity distributions, e.g., in gas turbulence [|l2l] or wide field scans with high 
spatial resolution [1 131 1 (in particular holographic image velocimetry displays diffraction limited 



spatial resolution 



1411 ). Both imaging applications, however, often require a complex instru- 



mentation, even if for 2D scans the less demanding speckle imaging velocimetry can be used lIlSll . 
Tracer experiments have also been reported on suspensions of interacting particles, but in general 
the choice of suitable systems is strongly restricted [16]. Laser Doppler velocimetry (LDV) on 
the other hand, which is a heterodyne light scattering technique, can be directly employed with 
pure samples and small particle sizes. LDV uses either real fringe or reference beam realizations 
linn . Using fringe techniques, two laser beams are crossed in the sample cell to produce a local 
illuminating grating. For sufficiently dilute suspensions the particle motion can then be traced 
via an oscillating scattered intensity of frequency proportional to the particle velocity. For large 
particle concentrations the analysis of the scattered intensity yields the heterodyne intensity au- 
tocorrelation function (lACF) or its Fourier transform, the power spectrum. Contributions of low 
frequency noise and the homodyne lACF (stemming from the self beat of photons scattered off the 
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sample) as well as the overlap of the two Fourier transform components of the heterodyne lACF 
[llSll (observed for low velocities) can be avoided by introducing a frequency shift between the two 
illuminating beams. Frequency shifting was first used for electrokinetic applications by Schatzel 
et al. to locally separate the diffusion and the small drift velocities in non-interacting particle sys- 
tems [19]. It may be noted that in that work also an alternative data analysis was employed to 
obtain the so-called amplitude weighted phase structure function from the phase of the scattered 
light [20] . The main general drawback of fringe techniques however is their strictly local velocity 
detection implying a tedious point-by-point determination of the actual flow pattern bill . 

By contrast in a reference beam set-up, a small portion of the illuminating beam is separated. 
It is re-combined at the detector place with the Doppler- shifted light scattered off a small vol- 
ume of the sample to act there as local oscillator and give rise to beats in the detected intensity 



2411 . The beat frequency is proportional to the Doppler frequency and therefore can be 



analyzed by considering either the lACF or its Fourier transform pair, the power spectrum. As the 
illuminating beam crosses the complete sample cell (integral measurement IISO), in principle one 
obtains information on the velocity distribution in a non-interacting sample, provided the scatter- 
ing particles are homogeneously distributed and all of same scattering power. Then the amplitude 
of the spectrum contains information on the number of scatterers moving with a given velocity, 
and for a known sample cell geometry the flow profile can be quantified. 

The situation of interacting samples, where a locally different structure will lead to different 
amplitudes has not been rigorously treated. The present paper will make a first step in that direc- 
tion. Further, in previous studies, the reference beam setup suffered from the restriction that the 
homodyne contribution centered around zero frequency could possibly overlap with the hetero- 
dyne scattering contribution containing the information on the particle velocities. In the present 
work, we use a straightforward extension of the approach, used by Schatzel et al. for fringe 
techniques, to integral reference beam techniques by introducing a frequency shift between illu- 
minating and reference beam. This shift is also accounted for by us in the calculation and analysis 
of the super-heterodyne signal. The general theory of dynamic light scattering from colloidal 
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25|]. Thehetero- 



particles is known since many years and has been summarized, e.g., in [|22L 
dyne (electrophoretic) light scattering technique for dilute systems of non-interacting particles has 
been discussed thoroughly in n24\ . The shifting method has been very shortly discussed in [|26ll . 
The present paper gives a theoretical derivation of the super-heterodyne intensity autocorrelation 
function and its Fourier-transform pair, the power spectrum. 
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Moreover, this derivation for the first time includes an account of the unavoidable polydisper- 
sity actually leading to a principal possibility to deal with differently structured systems. To be 
specific, the theoretical analysis of the super-heterodyne lACF reveals that the shape of the mea- 
sured power spectrum is usually connected to a convolution of the particle velocity distribution 
and the micro structure of particles, since the heterodyne part of the lACF contains an amplitude 
factor related to the local 'static' structure factor. Thus, for suspensions of correlated particles, 
the assumption of a homogeneous distribution of scattering power may be invalidated if the mi- 
crostructure of the system changes locally. It is a main point of the present paper to experimentally 
realize the solution to this problem as suggested by theory. Exploiting the dominance of incoher- 
ently scattered light at low scattering angles and strong particle interactions, we show that the 
spectra obtained from electro-osmotically sheared suspensions of fluid-like order still can be un- 
equivocally interpreted as velocity distributions. 

In addition we also observe an increase of the effective diffusion coefficient and the integrated 
spectral power with the electric field for the case of suspensions of interacting particles. The 
field-dependence of both quantities is approximately linear in the case of fluid-ordered systems 
under stationary flow, but more complex in the case colloidal crystals. Interestingly, the integrated 
power stays constant during shear-induced melting at constant applied field, suggesting that the 
additional spectral power is not induced by the structural change. 

While the physics behind these two effects are not yet fully understood, our tests clearly show 
that we have a versatile tool for quantifying the flow behavior in the fluid-ordered state, and for 
qualitatively complementing structural measurements in the crystalline state and in multi-phase 
flow. Given this, our approach should have several interesting applications. First, it could be 
quite useful to characterize deviations from simple Hagen-Poiseuille flow in non-Newtonian flu- 
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290. From a comparison of experimental results, it could 



be useful to separate the effects of shear, absolute drift velocity against the solvent and external 
fields on the suspension dynamics, thus enabling an extension of presently available theoretical ap- 
proaches, e.g., to investigations on the dynamic structure factor for interacting suspensions under 
flow [23]. Our method should be also well suited to study laning effects occurring in mixtures with 
different particle speeds [30] . And finally, it has already turned out to be useful for quantitative 



electrophoresis experiments in interacting suspensions yi 
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In what follows, we shall first explain the experimental set-up for integral super-heterodyne 
light scattering in the reference beam mode, then give a thorough derivation of the lACF and the 
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power spectrum for moderately polydisperse and strongly interacting charged- sphere suspensions. 
We then show that the assumptions made in our theoretical derivations can be well met in a typical 
experiment and verify the theoretical results with experiments on different experimental situations. 
We finally discuss the performance and scope of the super-heterodyne velocimetry method and 
give an outlook on open theoretical and experimental issues. 



II. EXPERIMENTAL 



A. Integral Super-heterodyne dynamic light scattering set-up 

Colloidal particle velocities have been obtained from a home-built super-heterodyne reference 
beam mode Laser Doppler velocimetry setup. Our instrument is a consequent extension of the 
heterodyne reference beam technique presented previously [llSl]. but in contrast to the latter con- 
tains a frequency shift between illuminating beam and reference beam. This ensures excellent 
discrimination of the electrophoretic signal from all other spectral components and low-frequency 
electronic noise. The set-up is shown schematically in Fig. [TJ The beam of a Nd-Yag laser (L, 
Aq = 532 nm) of circular frequency loq ~ 10^^ Hz is split into two beams, a reference beam (r) 
and an illuminating beam (i), using a beam splitter. Alternatively also a sinusoidal transmission 
grid (G) can be used (Sine Patterns, Pittsford NY). The two beams are made parallel by a lens 
Li and each beam passes through a Bragg cell (BC) and is frequency shifted by the circular fre- 
quencies LUi = 1.002 MHz and cUr = 1 MHz, respectively, with a positive relative Bragg shift, 
Acub = uji — ujr > 0, in the circular frequency corresponding to A/b = Ac<Js/27r = (2000 ± 5) 
Hz. The larger frequency is realized for the illumination beam. The beams are redirected into the 
sample cell (S) by lens L2 where they cross in suspension under an angle of 6 = 18.82. 

In addition to a polarizer (P) the detector side optics include a a set of apertures (A) to define 
the observed sample volume and an optical fibre to define the scattering vector q. The scattering 
vector is 

q = ki - (1) 

where kj and k^ are the wave vectors of the illuminating and scattered, respectively (c.f.. Fig. [l]). 
The modulus of the scattering vector is given as g = (AtcUs/Xq) sin(6/2) = 5.0 fim^^, where Ug 
the index of refraction of the solvent. For the present sign convention, q is proportional to the mo- 
mentum transfer from the photon to the scattering particle. For so-called integral measurements, 
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FIG. 1: Super-heterodyne velocimetry set-up. Top: Sketch of the rectangular flow-through cell and the 
used coordinate system. The plus and minus signs indicate the electrodes used to induce the combined 
electrophoretic-electro-osmotic flow. Middle: Top view of the super-heterodyne experiment (for a detailed 
description see text). Bottom right: definition of the scattering wave vector q. 

light scattered off the illuminated part of the sample is received within the acceptance angle of a 
grin lens. This lens is mounted at the inlet of an optical fibre leading to a photomultiplier used to 
record the intensity (D). Note that for the present definition of q, the light scattered by a particle j 
moving with a velocity vj is Doppler-shifted by the circular frequency lub = — q ■ Vj{x, y), where 
the velocity could be a function of the position. The frequency shift cud is positive, if the particle 
velocity shares an obtuse angle with the scattering vector q = g z, where z is the unit vector 
in z-direction. For example, in Fig. [T] a positive Doppler shift corresponds to particles moving 
towards the detector. Also the reference beam is directed into the fibre and superimposed with 
the scattered light. It therefore acts as a local oscillator, and gives rise to beats in the observed 
intensity which are analyzed by a Fast Fourier Transform analyzer (Ono Sokki DS2000, Com- 
pumess, Germany) to yield the power spectrum as a function of frequency / = tu/2n. As we will 
discuss, the spectrum is composed of a strong peak at zero frequency, a homodyne term stemming 
from the self -beating of the scattered light, which is localized around zero frequency, and a super- 
heterodyne part, which is Doppler-shifted relative to the Bragg frequency shift A/^. The latter 
term contains all the information on particle velocity and direction and is further evaluated. 

In our measurements, we have opted for a parabolic solvent flow profile. Such a profile can 
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FIG. 2: Horizontal x — z midplane: laboratory-frame velocity profiles v{x) = v{x,y = 0) and us{x) = 
us{x, y = 0) of particles and solvent, respectively. Vertical plane: corresponding probability density, P{v), 
of finding a particle of given velocity v. In our setup, the external field E, Ugo and q are aligned with the 
positive z-axis, whereas Vg points into the negative z-direction. 



be easily realized in a cylindrical tube, through which a suspension flows under a hydrostatic 
pressure difference. For strongly interacting suspensions, however, this approach often suffers 
from an insufficient stationarity of the flow, and an only moderate reproducibility [29]. Therefore, 
in the experiments discussed here, we have applied an electric field to move the particles and shear 
the suspension. Application of an external static electric field, E, pointing into the positive 2;-axis 
in a closed cell geometry (see Fig. [T]) allows for the development of an electro-osmotic flow profile 
of the incompressible solvent. Surface-released cations in the thin double layer at the negatively 
charged cell walls move towards the negative electrode with a velocity u^o > pointing along the 
positive z-axis (see Fig. [I]) and drag some solvent along with them. A central solvent backflow 
in direction of the negative 2-axis assures volume conservation of the incompressible fluid (Fig. 
121). For DC electric fields a stationary parabolic flow profile results in both cylindrical geometry 
and the mid-plane of a rectangular cell Oal . For optical reasons we here employ a rectangular cell 
as shown with its electrodes in Fig. [U The electro-osmotic flow profile, us{x), of the solvent is 
sketched in Fig. [2] for the x — z mid-plane located mid-cell of the rectangular cell at y = 0. Note 
that us{x) is equal to Ueo at the side walls, and due to the rectangular cell shape, the integral of 
us{x) with respect to x from one side wall to the other is zero due to volume conservation of the 
incompressible solvent. 

Superimposed on this solvent flow profile is the electrophoretic motion of charged colloidal 



8 



particles relative to the solvent with the electrophoretic velocity = He E, where /ig is the 
electrophoretic mobility. The constant electrophoretic velocity f g < of the negatively charged 
colloids points in negative z direction towards the positive electrode. The resulting flow profile, 



V{x) = Ve + Us{x) 



(2) 



of colloidal particles in the laboratory frame is thus equal to the solvent profile shifted by the 
constant amount \ve\ in negative z direction. Also drawn schematically in the Fig. [2lis the particle 
velocity distribution, P{v), where P{v) dv is the fraction of colloids with velocity in a infinitesimal 
interval dv around v. On assuming a homogeneous distribution of colloidal particles that remains 
unperturbed by the electro-osmotic flow, P{v) oc dx/dv [|l8n . 

Clearly the suspension is sheared at locally different rates by the underlying solvent flow. All 
colloidal particles and their electric double layers are subjected to the homogeneous electric field 
that couples differently to colloids and microions. If instead shear motion is introduced by a 
hydrostatic pressure difference and not by an electric field, the microions and the colloids are 
affected in the same way by the flow. If sedimentation is addressed in place of electrophoresis, 
gravity couples to the colloidal particles but not to the microions in the double layers and the 
solvent flow in sedimentation is not parabolic. Pipe flows are just one example of many possible 
shear profiles. For interacting colloidal particles, shear flow may lead to local phase transitions 
and, in concert with hydrodynamic interactions (HI), to a non-Newtonian flow behavior. This has 
been observed for both pressure and electro-kinetically driven flows [4, 5, ^7, 8, 21, 28.l29l. l3il 
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3411 . A fascinating richness of structural effects was there observed ranging from (partial) 



shear melting, to shear banding and shear induced structure formation, which all coupled back 
to the flow properties. In practically all cases deviations from the parabolic flow profile were 
observed. For the present paper we concentrate on one-phase systems of more or less pronounced 
fluid-like order and only shortly strive suspension flows with phase transitions. Consequently the 
theoretical section UlI] is concerned with the derivation of the lACF and the corresponding super- 
heterodyne spectrum for homogeneous systems in a stationary state. 



B. Samples and sample preparation 



One species of colloidal particles employed in this work are silica spheres synthesized by a 
modified Stober synthesis, showing a diameter of 2a = 251 nm and a relative standard deviation 
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in the size distribution of s = 0.08 determined from TEM-measurements. A representative TEM 
image is shown in the insert of Fig. [3l Suspensions of Si251 particles could be used in our 
measurements only in the fluid-like phase. Working at larger particle concentrations was inhibited 
by the onset of strong multiple scattering. The second species of spherical particles we have 
made use of is an industrial sample of Poly-n-Butylacrylamide-Polystyrene copolymer and a kind 
gift of BASF, Ludwigshafen. The PnBAPS68 particles have a diameter of 2a = 68 nm and 
a standard deviation of s = 0.05 (determined from ultracentrifugation) and carry an effective 
charge of Z = 420 obtained from conductance measurements. The particle size distribution 
deduced from ultracentrifugation is shown in Fig. [3l Aqueous suspensions of this species can be 
investigated up to large particle concentrations without too much multiple scattering, so that we 
have studied them also in the crystalline state. 

Both particle systems were prepared from pre-cleaned stock suspensions of large particle con- 
centration stored over a mixed-bed ion exchanger. By precise dilution with milli-q-grade water, 
and by determining the particle density, n, using static light scattering, the desired number den- 
sities have been adjusted and the resulting suspensions filled in a closed Teflon-tubing circuit 
connecting an ion exchange chamber, a reservoir under an inert Ar atmosphere (to add further 
suspension, salt or solvent), and a set of different measurement cells. These comprised a con- 
ductivity experiment, a cell for static light scattering, and the flow-throu gh c ell for the velocity 
measurements. Details of the used preparation circle are given elsewhere 11371 13 811 . We here only 
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FIG. 3: Particle size distribution from ultracentrifugation for PnBAPS68. The insert shows a TEM picture 
of Si251. The scale bar is 1 fim. 
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note that residual impurity concentrations on the level of the self -dissociation product of water are 
reached within a few hours. For measurements at elevated salt concentrations, NaCl solution was 
added after complete deionization. Continued cycling under bypassing the ion exchange column 
ensured homogenization, on controlling the salt concentration by conductivity measurements [|39D 
and the particle concentration by static light scattering with a residual uncertainty of about one 
percent. This preparation method is thus very fast and assures good reproducibility of adjusted 
suspension parameters. Furthermore, the preparation circuit permits to conduct several experi- 
ments simultaneously. This is of great importance to have a continuous control on the constancy 
of interaction parameters. 



C. Electrokinetic flow experiment 

Electrophoretic-electro-osmotic shear flow experiments were performed in closed quartz cells 
of rectangular cross section with depth x height = dx h = 1 mm x 10 mm (Rank Bros. Bottisham 
Cambridge, UK). As sketched in Fig. [Hon the upper left, the origin of our coordinate system is at 
the cell center with x \\ d,y \\ h and z \\ I. The optical part of the cell length, /, spans some 40 
mm and the effective platinum electrode distance is / = 85.5 mm as determined from calibration 
with an electrolyte dilution series. During the measurements the electrode chambers are sealed 
against the remaining preparation circuit by electromagnetic valves closed upon stop of pumping. 
To avoid accumulation of particles at the electrodes alternating square-wave fields of strength up 
to Emax = U/l = 80 V/cm were applied to fluid ordered systems. To ensure fully developed 
stationary flows, field switching frequencies /ac = (0.05 — 0.2) Hz were used. Measurement 
intervals were restricted to one field direction, starting 2s after field reversal and ending Is before 
the next field reversal. 



ni. THEORY 



In this section we give the essentials of the theory of super-heterodyne electrophoretic light 
scattering from suspensions of charge- stabilized colloidal spheres subject to a constant external 
electric field. We proceed from non-interacting particles to low-salinity suspensions of strongly 
interacting particles. To our knowledge, the super-heterodyne electrophoretic scattering theory for 
strongly interacting colloidal particles has not been presented so far. 
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For a stationary and spatially homogeneous system, we derive the general expression for the 
measured intensity autocorrelation function (lACF) and its related power spectrum (PS) in terms 
of (partial) dynamic structure factors. The structure factors characterize the density fluctuation 
correlations of colloidal spheres in the presence of an external electric field. It will be discussed 
on how the electrophoretic velocity of the spheres may be deduced from small wave number 
measurements on weakly poly disperse systems. 



A, Super-heterodyne scattering theory 

The super- heterodyne scattering setup has been shown schematically in Fig. [Tl We assume in 
the following that the colloidal spheres are optically isotropic, and that the Rayleigh-Gans-Debye 
description of weakly scattering particles applies where the electric fields scattered at each point 
inside a sphere can be linearly superimposed. The following general discussion does not include 
electro-osmotic flow effects which are of importance in closed cell geometries. The influence of a 
non-zero macroscopic solvent flow profile on the measured super-heterodyne power spectrum will 
be discussed in subsection ini F[ 

Let the Bragg-shifted illuminating electric field be a plane wave of the form 

Ei(r, r) = fii E° exp{i (k^ ■ r - [ujq + Ui] r)} , (3) 

where fij is the unit polarization vector of the incident field. Under the assumptions made above, 
and for the standard V-V geometry of polarized scattering, the field strength, -Es(q, r), of single- 
scattere d lig ht from spheres in the scattering volume follows from standard scattering theory as 
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m 



N 



Es{<l,r) = exp{-i{uJo + uJi)T}c^fi{q) exp{i(i- riir)} , (4) 

1=1 

with q as defined in Eq. ©• Here, fi{q), is the scattering amplitude of a sphere labelled by /, 
whose center at time r is at the position ri^r). For a homogeneous sphere of radius a;, volume V/ 
and refractive index vi, the scattering amplitude is given by //(g) = 2(z/; — z/^) bi(q)/us, with 
the form amplitude bi{q) = 3ji{qai)/{qai), where ji is the spherical Bessel function of first order. 
The apparative constant 

_ exp{i\K\R} 
' - A^i^ 

depends on A, kg, the large sample-to-detector distance -R, and the field amplitude, E^, of the 



illuminating light 11221 l23h. Its explicit form is of no concern, constant throughout the experiments 
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and thus will not be displayed any more. The corresponding electric field strength of the reference 
beam at the detector is 

Er{T) = exp{~i {loo + Ur) t} . (6) 

The reference field is independent of the scattering wave vector and the particle configuration. 

In the (super-)heterodyne setup, the intensity autocorrelation function of the mixed electric 
fields is analyzed by a real-time auto-correlator. For a stationary system, the super-heterodyne 
mixed- field intensity autocorrelation function (lACF) is given by 

ahet(q,r) = {\E,iq,T) + EriT)f-\E,{q,0) + Er{0)f), (7) 

where r is the correlation time, which in our experiments is 10 — 100 ms or larger. The angular 
brackets denote a time average or likewise, for a stationary and ergodic system, an ensemble 
average. The lACF depends on the orientation of q, since spatial isotropy is broken by the applied 
electric field. Carrying out the multiplications gives sixteen terms, ten of which vanish due to 
the assumed spatial homogeneity, and three of which are time-independent. The remaining three 
terms determine the time dependence of the lACF. To see this, we rewrite Eq. ([7]) as 

ahet(q, r) = ([/,(q, r) + /,(r) + Er{r)E:{q, r) + i?;(r)E,(q, r)] x 

[/,(q,0) + 1,(0) + E,(0)E:(q,0) + E;(0)E,(q,0)]) , (8) 

with /s(q, r) = i?*(q, r) ■ Es{q,T) and /r(T) = E*{t) ■ Erij). The asterisk denotes the oper- 
ation of complex conjugation. The reference beam quantities, and E^, are independent of the 
ensemble average over the particle configurations. Moreover, the averages over products of an 
odd number of factors Es{q) and/or El{q) vanish for g > 0, since spatial homogeneity requires 
the wave vector sum to be zero. For the same reason, (-E's(q, r)Es{q, 0)) and {E*(q, r)_E'*(q, 0)) 
are identically zero. Using these properties of spatial homogeneity, we multiply the first element, 
/s(q, r), in the upper bracket of Eq. ([8]) with the elements of the lower bracket, keeping terms 
only which are non-zero after averaging. On repeating this multiplication for the remaining three 
elements in the upper bracket, we obtain the intermediate result 

Cshet(q,r) = Ir[Ir + 2{Isiq))]+2Rc{E:.{T)EriO) {Es{q,T)-E:{q,0))} + 

{\E,iq,T)f-\E,{q,0)f) , (9) 

showing that three out of six non-zero terms are time-independent, and add up to a constant base- 
line Ir [/r + 2 (/^(q))], where = \E^f and /^(q) = \Es{q,0)f. To make the frequency depen- 
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dence of Cshet(q, t) explicit, we redefine Es{q, r) now as 

N 

EsicL, r) = cJ2 fM) exp{^q ■ r^r)} . (10) 

1=1 

With this redefinition, we obtain 

Cshet(q,r) = Ir [Ir + 2 (/.(q))] + 2/,Re [^^^(q, r) exp{-2 Acj^r}] + ^/(q, r) , (11) 

for the super-heterodyne scattering function, where Acj^ = uji—ujr > ^- The standard expression 
for heterodyne scattering is recovered for zero Bragg shift [25]. Note that Cshet(q)T) does not 
depend on ujq since, due to the mixing of -Es(q, r) and E*{t), it only contains the low-frequency 
beats. 

Information on the particle motion is embedded in the super-heterodyne part of the mixed-field 
lACF, given by the electric field- autocorrelation function (EACF) of the scattered light 

(7s(q,r) = (E,(q,r).i5;:(q,0)) . (12) 

The mixed-field lACF contains also a homodyne part, equal to the homodyne intensity autocorre- 
lation function, 

gi{c^,T) = Chom(q,r) = (|i?.(q,r)|2 ■ |i?.(q,0)|'> (13) 

which is of order 0{E^). We further note that 5's(q, 0) = (/s(q)), and 5'/(q, 0) = (/^(q)). In 
standard heterodyne experiments, (J^) ^ J^, and the homodyne contribution is negligibly small. 
Eq. (fTTI) reduces then to 

Cshet(q, r) = + 2Ir (/.(q)) ( 1 + Re [ M'i, r) exp{-tAujBr}] ) , (14) 

with the normalized EACF of the scattered light, 

gEiq, r) = , (15) 

defined such that ^£;(q, 0) = 1. The EACF is complex- valued for a non-isotropic system in an 
external field. 



B. Scattered field with Gaussian statistics 



So far, we have not assumed that the scattered field is a complex central Gaussian random 
variable, an assumption frequently made in dynamic light scattering theory. This assumption is 
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certainly valid for dilute suspensions with weak particle correlations, where it becomes a conse- 
quence of the central limit theorem. It is violated, however, in very dense suspensions of glass-like 



or gel-like character that show dynamic heterogeneity [411]. and in live sperm cells . For a scattered 



field of Gaussian character, Eq. (fTT)) can be derived alternatively using Wick's theorem for the av 



erage of products of central Gaussian random variables [|23l] . Another consequence of Wick's 
theorem is that for a Gaussian scattered field, 5'/(q, r) includes no information on the suspension 
dynamics, not contained already in (7s(q, r), and that they are related by the Siegert relation 

?/(q,^) = 1 + \9E{q,T)f , (16) 

with ^/(q, 0) = 2. The mixed- field lACF for a scattered field of Gaussian statistics is thus 

Cshet(q,^) = [/r + (/.(q))]^ + 2J^(J^(q))Re[^£;(q,r)exp{-iAu;Br}] 

+ (/s(q))'fe(q,r)|' (17) 

with the normalized scattered-field lACF, 

9i{(l,r) = (18) 
(^4q)) 

defined such that ^7 (q,0) = (/^(q)) /(/,(q))^ > 1. For a non-Gaussian scattered light field, 
however, one has to deal with two genuinely different correlation functions ^'^(q, r) and gi{q, r). 



C. Super-heterodyne power spectrum 

When a spectrum analyzer is used in place of a time correlator to analyze the superposition of 
scattered and reference light, the power spectrum is determined in place of the mixed-field lACF. 
The power spectrum is the temporal Fourier transform of the mixed-field lACF, namely 

2 /"OO 2 /'OO 

Cshet(q, cu) = TT dr exp{iujT}Cshet{(l,r) = - rfr cos{u;r}Cshet(q, , (19) 

where the second equality follows from the stationarity property Cshet(q, t) = Cshet(q, — t). Sta- 
tionarity implies further that (7s(q, r) = gEi<l, — t)* and (7/(q, r) = 5'/(q, — r). The mixed-field 
lACF and its power spectrum are measurable quantities and therefore real- valued. The power spec- 
trum is symmetric in cu, that is it consists of two sub-spectra extending symmetrically to positive 
and negative frequencies. Furthermore, it fulfills the sum rule 

/oo 
dujaU^,u) = ahet(q,r = 0) = + 4/,(/,(q)) + (/.(q)^) . (20) 
-oo 
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The 'static' quantity (/s(q)) = (/s(— q)) in a driven, non-equilibrium system still depends on the 



hydrodynamic interactions ll42ll . 



D. Non-interacting colloidal spheres 

Consider a very dilute charge-stabilized suspension, which includes a sufficient amount of 
added salt ions, so that the colloidal spheres are non-interacting. Suppose that the spheres 
in the scattering volume form m components, with Na particles in component a E {1, ...,m}. 
The position vectors of non-interacting spheres are independent so that 

m 

(C(q)) = |cr5^iV„/2(g), (21) 

a=l 

and 

&°(q.r) = 2S^i?^M^. (22) 

Here, Xa = N^/N is the molar fraction, fa{q) the scattering amplitude, and 

G„°(q,r) = (exp{zq- (r^r) -r^O)}) , (23) 

the self-intermediate scattering function of component a. The label (°) indicates properties of a 
system of non-interacting particles. 

The self-intermediate scattering function is the spatial Fourier transform of the fundamental 
solution Ps{v, t\yq), with Ps(r, 0|ro) = 5{y — fq), of the Smoluchowski equation, 

§^Ps{r,r) = DoV'Ps{r,T)-^r:■VPs{r,r), (24) 

which is used to describe a single microion-dressed charged colloidal sphere of diffusion coeffi- 
cient Dq in an unbounded and quiescent fluid subject to a constant electric field E. The colloidal 
macroion responds to this forcing by acquiring a terminal mean drift velocity relative to the quies- 
cent solvent 

v: = /i:E, (25) 

once the polarization of the microionic electric double layers (EDL) in response to the external 
field has acquired a steady state after the relaxation time tdl = (/^^-Dei) \ which is typically 



10 s for a 1 mM electrolyte solution [|43ll . Here, /x° is the electrophoretic mobility of a colloidal 



16 



macroion in isolation, and and are the Debye screening length and mean diffusion coeffi- 
cient, respectively, of electrolyte ions. The assumption underlying the use of the forced diffusion 
Eq. (|24)) is that the particle diffusion and the electrokinetic effects on the EDL caused by the ap- 
plied field are independent. This assumption is reasonable for a weak field which can be treated 
as a first-order perturbation (linear electrophoresis). The field-independent mobility of a spheri- 



cal colloid can be calculated using standard electrokinetic mean-field theory (see, e.g., 11431 
or extensions which account for microion correlations [45]. These theories provide a relation 
between the experimentally determined electrophoretic mobility and the zeta potential or related 
electrophoretic colloid charge. 

The fundamental solution of the single-colloid Smoluchowski equation is 

P,{r, T > 0|ro) = -—^-jH { - ^^-^^^r^l ' (26) 
{An Dory ^ ^DqT J 

stating that r — fq — v°r is a central Gaussian random variable so that (7/(q, r) is related to 5'^(q, r) 

by the Siegert relation. 

Eq. (|26l) applies also to a polydisperse system of statistically independent spheres of varying 
electrophoretic drift velocities. Electrophoresis differs from sedimentation in that the total force 
exerted on the suspension is zero due to overall electro-neutrality. This has consequences for the 



range of the HI, which is shorter in the case of electrophoresis [|46L I47|l . 

The self-intermediate scattering function of an isolated sphere of diffusion coefficient Dq = 
ksT / (Gnrjoo), which drifts with the mean velocity v° = /i" E, follows from Fourier-transforming 
the fundamental solution, 

G°(q,r>0) = j d{r -Tq) exp{iq- {r -ro)}P,{r,T\ro) 

= exp{2q- v°r}exp{-g2L>or} , (27) 

where the sign convention in the Fourier transform has been made consistent with the one in Eq. 

In summary, the normalized FACE and lACF of a monodisperse suspension of non-interacting 
colloidal spheres with equal electrophoretic drift velocity v° is given, for all r, by 

giiq,T) = exp{^q■^r:T}exp{-q^Do\T\} (28) 
?/(q,r) = l+exp{-2q'Do\r\}, (29) 

which shows that the homodyne part of the mixed-field lACF is independent of v° as long as no 
differential motion between two spheres due to different drift velocities occurs. Polydispersity in 
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the particle mobility would introduce such a differential motion between different components a 
and P and thus be detectable as an oscillation in ^/(q, r) ll48n . Continuing with the case of equal 
drift velocities and using Eq. (fTTI) . the mixed-field lACF becomes 

<^shet(q, r) = [Ir + (/:(q))]' + 2Ir (//(q)) COS [(A^B - q ■ v:) r] expj-g^Do |r|} 

+ (//(q))'exp{-2g2Z}o|r|}, (30) 

where //(q) = \c\'^N p{q). Suppose now that {If) < /,.. Then, 

Cshet(q, r) = I? + 2/, (4«(q)) + 2/, (//(q)) cos [(Aa;^ - q • v°) r] expl-g^D^ |r|} , (31) 

i.e., provided the homodyne part is negligibly, Cg^etlq, t) is an exponentially damped cosine in r, 
shifted upward by a time-independent baseline, and characterized by a 

1 



decay time: r, 



q^Do 



period: T = - — -. (32) 

\Aljb - q ■ v°| 



Thus, a measurement of the mixed-field lACF provides a simultaneous determination of the dif- 
fusion coefficient and the electrophoretic mobility, given in our setup by = |(q ■ v°)|/(g E). 
Super-heterodyning using Bragg cells has two effects not present in a conventional heterodyne 
setup: First, through the appearance of a non-zero Alub, the oscillation period T can be short- 
ened relative to the decay time r^, which allows for an improved resolution in v^, and second, the 
method is sensitive to the sign of v°. A drift velocity parallel to q || E || z or, more generally, a 
velocity sharing an acute angle with q, enlarges the period of oscillations, whereas antiparallel 
motion causes the period to decrease. In our setup, the negatively charged colloids move antipar- 
allel to the q-axis (i.e., — q ■ Vg = g |t>g | > 0), so that a reduction in the period T is achieved. The 
resolution can be further improved by making larger, that is by choosing the scattering angle 
as small as possible. For typical colloids, r^^ = Au ~ 10"^ Hz, so that the symmetric frequency 
broadening due to the isotropic diffusive motion of spheres is comparable to the frequency shift 
by the Bragg cells. 

Using Eq. (fT9l) . the mixed-field power spectrum is obtained as 



+ 



TT 



+ 



CO + [AcuB - q • v^])' + (q^Dof (to - [Acub - q ■ v°])' + (q^ Dof 

\2 o^2 



(«3)>:^^. (33) 
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It is the sum of an irrelevant singular term at u = due to the non- subtracted baseline contri- 
bution, and two symmetrically shifted Lorentzian curves of half-width at half-height Dq, cen- 
tered at = ± [Acub — q • which is the interesting super-heterodyne part, and an unshifted 
Lorentzian of double- sized half- width at half-height 2q'^Do, centered at = 0. The last contribu- 
tion is due to the self-beating of the scattered light (homodyne part) and contains no information on 
the electrophoretic motion. The two heterodyne Lorentzians are Bragg-shifted away from the ori- 
gin by Alub, so that they can be better distinguished from the central homodyne Lorentzian. Since 
q is parallel to E = £'z in our setup, and since the negatively charged colloids move antiparallel to 
the field, the centers of the super-heterodyne Lorentzians are separated from the frequency origin 
by Acub + q fx^E, so that this separation increases with increasing field strength. The frequency 
resolution improves with decreasing q and Dq, and increasing field strength E. Note that the fre- 
quency shift in the mixed-field power spectrum corresponds to the sinusoidal modulation in the 
mixed-field BACK 

When the electrophoretic mobilities vary from particle to particle possibly due to some charge 
polydispersity, the super-heterodyne part of the mixed-field lACF in Eq. (fTT)) is generalized to 



^shet(q, ^)lpart = 2/. (//(q)) Rc ( ?|(q, t) ) exp{-z Ao^^t} 



(34) 



where {■ ■ denotes an average over the mobility distribution. A narrow distribution character- 
ized by the mean mobility /Zg and the relative standard deviation o^^ , can be approximated by a 
Gaussian, yielding [|48ll 

^^shet(q,^)lp^. = 2/.(//(q))exp{-g2 



Do|r| + ^(7Ie^r)Vj^ 



}cos[{AuB + qJleE)T],05) 



for non-interacting colloidal spheres. Hence the super-heterodyne part of the lACF decays faster 
than exponentially with a rate that increases with increasing mobility polydispersity cr^. The cor- 
responding power spectrum has a shape similar to a Lorentzian centered at Aujb + qJ^eE, but is 
broadened symmetrically by an amount that increases with increasing a^^ . 

A remark on the sign convention for q made in Eq. ([T]) is in order here. In fact, we could 
have equally well used the definition q = — (k,; — kg) for the scattering wave vector, which 
amounts to replacing q by — q in Eq. dH). The only consequence is that (7£;(q, r) must be replaced 
by q, r) = gEi<l, t)* in all the following expressions, on noting further that (?/(— q, r) = 
5f/(q, r) and q)) = (/^(q)). For non-interacting particles, this leads to the replacement of 
Alub — q ■ Vg by Alub + q ■ Vg in Eqs. ([30] -[33]). This makes it obvious that the outcome of the 
experiment is independent on the sign convention used for the scattering vector. 
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E. Interacting colloidal spheres 

Information on the colloid velocities is embodied only in the super-heterodyne part of 
Cshet(q, t) that involves gEi<i, t)- The normalized EACF of light scattered from interacting 
spheres is given by 

^ , , EZ=iM^)fM (exp{^q. [r,(r) - r,(0)]}) 

5fE(q, T) = . (36) 

Ei,i=i fiil) fAl) (exp{iq ■ [r; - r^]}) 

Consider first an ideally monodisperse suspension of identical spheres. Then, 

c/i? q,^ = ^, . , (37) 
S{q) 

and 

(/,(q)) = \cfNf{q)S{c[). (38) 

Here, 

^(^'^) = ^ exp{zq- [r,(r) -r,(0)]}|) (39) 

is the complex-valued, steady-state dynamic structure factor. The corresponding steady-state 
structure factor is S'(q) = ^(q, 0). 

For a m component mixture, it follows from Eq. (l36l) that 

^, . E^/3=i (a;aa;^)^^Va(g)//3(g)'^a/3(q,r) 

gE{<l,T) = , (40) 

with the partial dynamic structure factor of a — P colloid pairs defined by 

^-M^' = ( E E ■ - } ) ' (41) 

where points to the center of a sphere / belonging to component a. The partial dynamic structure 
factor can be written as the sum of a self- and distinct part, 

Safi (q, r) = 6o^f3 (q, r ) + S^^ (q, r) , (42) 

where the self-intermediate scattering function of an a-type sphere is defined such that Ga(q, 0) = 
1. On invoking the definition 

^ m 

Sm{(1,t) = === V (XaXf^)^^^ f^{q)f,3{q)Sai3{q,T) (43) 
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of the measurable dynamic structure factor in a polydisperse system, and 

m 

Mq) = Y^xJMf. (44) 



a=l 



for any integer p, Eq. (|40l) is rewritten in more compact notation as 

withS'Af(q) = S'A/(q,0). 

Suppose now that the m components differ only in their scattering amplitudes {/a(q')}, but are 
otherwise identical. For such an interaction-monodisperse system, the partial dynamic structure 
factors reduce to [|2|] 



SaM^ r) = 5,^3 G(q, r) + {x^xpf^ [^(q, r) - G(q, r)] , (46) 

where (^(q, r) and ^(q, r) are, respectively, the self-intermediate scattering function and the dy- 
namic structure factor of an ideally monodisperse system. Therefore, for pure scattering polydis- 
persity, the EACF of the scattered light reduces to 

QEici^r) = |c|'iV7^(g)5^(q,r), (47) 

with the so-called decoupling (approximation) dynamic structure factor 

5z5(q, r) = X{q) ^(q, r) + [1 - X{q)] 5(q, r) , (48) 

as an approximation to S'a/ (q, r), and the decoupling scattering factor 

X{q) = 1 - £ilL . (49) 
P{q) 

For a suspension of slightly size-polydisperse homogeneous spheres made of the same scattering 
material, X{q) simplifies to 

X(g) ^ 9s^ (50) 



provided that s < 0.1 and qa < 1.5, where cr = 2a is the mean diameter. The polystyrene sphere 
samples explored in this work are only slightly size-polydisperse, characterized by the relative 
standard deviation of s ~ 0.05. For these systems, the assumption of interaction-monodispersity 
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is thus quite reasonable, and it implies for 9 ^ 1, and in the experimentally probed small-g 
range, that 

QeM ^ \cfNf\q) [9s2G(q,r) + 5(q,r)] , (51) 

where terms of 0{s^) and 0{s^ x S'(q)) have been ignored. Here, /(g) is the scattering amplitude 
of ideally monodisperse spheres of diameter a. 

For values of q well below the position, qm, of the principal peak, the zero-field S{q) of deion- 
ized suspensions is very small in comparison to one, since the osmotic compressibility is very low. 
In these deionized fluid systems, qm ~ 27m^/^ = (67r^0)^^^ /a, where is the volume fraction of 



spheres ll49ll 



On assuming that S'(q) remains very small for |q| -C qm, also in the presence of a weak external 
field, as it can be expected in the linear response regime, we obtain approximately 

^ij(q,r) ^ G(q,r) , (52) 

provided that q <^ q-m and S'(q) ^ 9 s^. The scattered-field EACF in the super- heterodyne part 
of Eq. ([TT]) is then determined essentially by the self-intermediate scattering function, without 
significant contributions from ^(q, r). 

We point out here that the field-free coherent scattering contribution, S{q, r), to Soiq, t) de- 
cays faster than G{q, r) at low q, even at longer times when the slowing influence of particle caging 
comes into play. This is due to the large value of the (field-free) collective diffusion coefficient, 
Dc = DoH^q — > 0)/S{q 0) in comparison to Dq, where H(q) is the so-called hydrodynamic 
function [49]. In fact, G{q <^ qm,T) ^ exp{— W^l^)} decays slowly at small q since the parti- 
cle mean-squared displacement, W{t), has a long-time slope equal to the long-time self-diffusion 
coefficient, that is significantly smaller than Dq. For an example, near the freezing transition, the 
long-time self-diffusion coefficient amounts to only about ten percent of Dq. Therefore, G(q, r) 
is expected to dominate 5'(q, r) even in presence of a weak external field. 

It is a difficult task to calculate the dynamic structure factor, ^(q, r), of interacting macroions 
in the presence of an electric field since this quantity depends on all ionic degrees of freedom. In 
fact, to our knowledge, no general theoretical or simulation results are available to date for the 
dynamic structure factor of interacting charged colloids in an external electric field. Results are 
known only for very simplifying models of driven colloid systems. To give an example, the sta- 
tionary dynamic structure factor of a homogeneous and ideally monodisperse suspension of slowly 
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sedimenting colloidal spheres where (quite unrealistically) hydrodynamic interactions have been 
ignored is given by ^(q, r) = S'eq(g, r) xexp{zq- VgedT}. Here, Vsed is the sedimentation velocity 
measured in the laboratory frame, which is corrected for the solvent backflow originating from the 
presence of the container bottom, and S'eq(g, r) is the equilibrium dynamic structure factor of the 
non-driven system. Many-body hydrodynamic interactions and polydispersity in a suspension of 
particles subject to a constant forcing cause a change of the equilibrium microstructure into a new 
non-equilibrium state so that the use of a simple Galilei transformation becomes invalid. Li elec- 
trophoretic experiments, the electro-kinetic coupling of ions and possibly existing electro-osmotic 
solvent flow are additional complications that add to the complexity in calculating ^(q, r). 

To rationalize heuristically how the electrophoretic mobility of strongly interacting particles 
may be determined from a low-g measurement of G'(q, r), we make the assumption, just like for 
non-interacting particles, that diffusion and electrophoretic drift are uncoupled in a weak external 
field. We impose this assumption on a coarse-grained level only, where distances larger than 
Ax ^ '2,11 / qm are resolved, corresponding to long correlation times r Tm = (g^ Do)~^. 
On ignoring in addition polydispersity effects, the coarse-grained single-particle mean current, 
jg(r, r), associated with coarse-grained single-particle density Ps(r, r), can be decomposed as 

J, (r , r ) = (0) Vp, (r , r ) + Ve (0) (r , r) . (53) 

In this phenomenological description, -Ds(0) is the unperturbed long-time self-diffusion coeffi- 
cient, and Ve(0) is identified with the long-time mean electrophoretic velocity of the spheres, hi 
conjunction with the continuity equation, 

^p,(r,r) + V-J,(r,r) =0, (54) 

the single-particle eq. (|24l) is recovered, but Ds{(f)) is replacing now Dq, and Ve(</)) replaces the 
single-sphere electrophoretic velocity v^. The self-intermediate scattering function is here given 
by 

Gih\ qm,r :^ Tm) = (Ps(q,i")Ps(-q7 0)) exp{?q- Ve(0)r} exp{-g^L's(0)r}. (55) 

For deionized aqueous suspensions of PnBAPS68 spheres the system parameters are cr = 68 nm, 
^ 10^^ T = 293.15 K, z/,(A,T = 20C) ^ 1.33 (water), A = 488 nm in the static structure 
factor measurements for zero external field, and A = 532 nm in the electrophoresis setup with 
|E| ^10- 120 V/cm. This gives Ta = / Dq = 0.18 ms, Dq = 6.3pmVs, {q*f Dq = 156 Hz 
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corresponding to r^* = 6.4 ms, and q* /q-m ~ 0.4 where q* = (AnUs/X) sin{0/2} = 4.97 ;um~^ 
is the wave number used in the electrophoretic setup. The correlation times are in the range of 
10 — 100 ms or larger so that long-time particle motion is probed. 

The assumption on the decoupling of diffusion and drift is flawed when the electric field is so 
strong that -Ds(0) and Ve(</') become field-dependent [50]. To test the validity of Eq. (|55l) in a cell 
geometry with negligible osmotic flow requires to verify whether the super-heterodyne peaks in 
Cshet(q5 are of a Lorentzian shape and, if this has been confirmed, to determine the diffusion 
coefficient from the measured half-width at half-height, Act;i/2, according to Ds = AuJi/2/q'^. This 
coefficient should be substantially smaller than Dq in the case of a deionized system, provided it 
can be identified with the unperturbed long-time self -diffusion coefficient. The short-time self- 
diffusion coefficient in these systems, on the other hand, is only slightly smaller than Dq, with the 
concentration dependence given by (l ~ ctt^''^^), where ^ 2 — 3 [|49ll . 

Finally, to check measured velocities, the electrophoretic mobility is determined by the center 
of the heterodyne Lorentzian contribution to the power spectrum. It can be compared to theoretical 
predictions obtained from different methods. Quite recently, substantial progress has been made in 
developing novel numerical schemes that allow to calculate the electrophoretic mobility of charged 
colloids at non-zero concentrations. These schemes are based, respectively, on a smoothed surface 
method, where the sharp colloid-fluid interface is replaced by a diffuse interface 115 ill , on a fluid 
particle dynamics (FPD) method with charge-densities included ^2\, and a hybrid simulation 
method that combines the Lattice Boltzmann fluid description with a Langevin equation treatment 



of the colloidal particles [|53L 



54|l . A semi-analytic method to compute the electrophoretic mobility. 



/ie(0), in dense systems based on a simplified mode-coupling scheme applied to the Primitive 
Model, has been developed very recently by one of the present authors [55]. All these methods 
allow, in principle, to explore the electrophoresis of charged colloids with strongly overlapping 
electric double layers. Also in analytical theory, considerable progress was reported in extending 
the standard electrokinetic model of dilute systems, and dense systems with weakly overlapping 
double layers [44, 5^ based on the single-macroion cell model, to more dense systems of particles 
with strong electrostatic interactions [i57J. 
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F. Solvent flow profile 

The scattering theory described so far assumes a uniform motion of all spheres. This is exper- 
imentally realized in cells where the electrodes have a parallel plate capacitor geometry and are 



placed far away from any container wall to avoid electro-osmotic solvent flow pSi] . Cations in 
solution tend to concentrate at the negatively charged (glass) walls and move towards the cathode, 
dragging solvent (water) with them. This flow establishes a hydrostatic pressure which tends to 
force the center of the fluid to flow towards the anode, establishing thus a parabolic back flow in 



the cell [1590. It is an open question whether the osmotic flow profile is affected significantly by 
a non-zero concentration of strongly interacting colloidal macro-ions (usually negatively charged) 
which move towards the anode and drag also solvent along with them. In a first approximation, 
however, this feedback of the colloids on the fluid flow may be ignored. Moreover, as discussed in 



[|60ll . to built up an electro-osmotic profile takes quite a large time of the order 0.1 — 1 s, and this 
time should be accounted for in experiments with an alternating electric field. In our setup, the 
switching period is 15 s or larger so that the mobility is measured in the steady-state with a fully 
developed flow profile. 

The electro-osmotic flow profile in the x — z midplane at ?/ = of a rectangular vessel, for 
the solvent streaming along the z-axis with the velocity us{x = ±d) = u^o very near the cell 
walls, has been calculated by Komagata [36|. On assuming that this profile is macroscopically 
unperturbed by the interacting colloids, and that the electrophoretic particle velocity, < 0, 
measured relative to the solvent remains constant under the shearing solvent motion, we obtain 
using Eq. (O, 



V{x) = Ve + Us{x) = Ve + Ue 



1-3 



(56) 



_2 - 384/(7r5K), 

for the velocity profile of the particles in the laboratory (cell) frame. Here, K = h/d, and ±d 
and ±h are the locations of the cell walls in the x and y-direction, respectively. We see here 
that the solvent profile us{x) is proportional to the electro-osmotic velocity, Ueo > 0, which for 
small fields grows linearly in E. It is thus characterized by the zero-field mobility, fi^o = Meo/|E|, 
of microionic cations in the mid-plane close to the cell walls. Again we have assumed equal 
mobilities and thus equal Ve for all particles. A distribution of velocities can conveniently be 
introduced at this place, when a coupling of differential velocities and structure is neglected. 

In our experiments the aspect K is equal to 10. For the theoretical interpretation of the exper- 
imental power spectra of interacting colloids, we assume now that the homogeneous distribution 
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of particles remains unperturbed by the locally shearing motion of the solvent flow profile so that 
the particle velocity distribution derives from the Komagata profile according to P{v) dv oc ndx, 
with n constant. This implies, in turn, that the solvent profile is macroscopically unperturbed by 
the interacting colloids. The explicit form of the normalized P{v) derived from the Komagata 
profile is given in jlSll and will not be repeated here. In fact, for a sheared suspension of (neu- 
tral) rigid colloidal spheres, Brownian motion will give rise to a uniform particle distribution. A 
flow-induced cross-stream migration is observed, however, for flexible particles like polymers and 



poly electrolytes. At 



distribution ||61 , 



62L 



larg e Peclet numbers, this migration can lead to a non-uniform center-of-mass 
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G. Spectral power for a non-quiescent solvent 

In making use of the above considerations, we now proceed to a prediction of the power spec- 
trum shape for colloidal particles in an electro-osmotic flow experiment in a closed cell. We follow 



jlSll and use the convolution integral 

Cstt(q, /) = r^^^" dvP{v) Cl,(q, /; v, Deff) , (57) 

as an approximation to the measured power spectrum C^^^{q, /). Here, C°jjgj(q, /) is of the form 
of the single-sphere power spectrum given in Eq. (l33l) . with / = tu/(2n). The integral over the 
normalized particle velocity distribution function P{v) derived from the Komagata solvent profile 
extends over a finite velocity range depicted in Fig. [21 The velocity-averaged power spectrum 
depends also on u^o through its dependence on P{v). Note that (super)-heterodyne scattering is 
sensitive to particle velocities taken relative to the laboratory (cell) frame. 

We summarize here the assumptions underlying Eq. (l57l) : a) the particle electrophoretic ve- 
locity is neither coupled to the diffusive motion nor to shear and the stationary parabolic flow 
profile is unperturbed by the particle interactions (if present); b) the mean particle density remains 
unperturbed by the shearing motion of the solvent and the applied electric field; c) interaction 
polydispersity and polydispersity in the electrophoretic velocities are disregarded, and d) density 
inhomogeneities are not considered irrespective of structure formation in strongly interacting par- 
ticle systems. 

Furthermore we assume e) on exploiting the dominance of incoherent scattering at low q, that 
the power spectrum of a homogeneous and quiescent suspension of strongly interacting colloidal 
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spheres can be described approximately by the single-particle power spectrum C^^^^^{q, f) in Eq. 
(l33l) . Here, we make two replacements. First we substitute Dq by Df.ff{(f)), where in the absence 
of coupling between diffusion and shear we expect that Deff((t> ) = Ds{(f)), while at higher shear 



rates larger values of Dgff may be expected [|42l . I65L l66L l67i l681]- Second, we replace by v{(f)). 
Likewise, the 'static' prefactor (/°(q)) multiplying the relevant super-heterodyne part in Eq. (133]) 
should be replaced by (/^(q)) oc np{q) [9 + 5'(q)] for a weakly size-polydisperse system. 

The power spectrum predicted by Eq. (l57l) . which is a continuous distribution of velocity 
shifted Lorentzians, is shown in Fig. |4l It has the typical super-heterodyne form consisting of 
two symmetrically displaced heterodyne wings at about ±A/b, where A/^ = AcoB/C^n), a 
sharp zero-frequency background term indicated by the two narrow vertical lines, and a central 
homodyne portion symmetrically distributed around zero frequency and unaffected by the velocity 
distribution. The form of the particle velocity distribution depicted in Fig. [2] is recognized easily 
in the shape of the super-heterodyne wings, broadened slightly due to diffusion of the particles. 
Experimentally, only the positive-frequency heterodyne wing of the symmetric power spectrum is 
resolved, which includes all the information on the velocity distribution and, thanks to the super- 
heterodyning, is well separated from the non-interesting homodyne part, and there is no overlap 
of the symmetrically displaced super-heterodyne wings of the power spectrum. 

For the analysis of the experimental power spectra in terms of the theoretical spectrum de- 
scribed by Eq. (|57l) . the latter is fitted to the experimental power spectra data (super-heterodyne 




FIG. 4: Theoretically expected power spectrum according to Eq. (l57l ). 
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wing) using four independent fitting parameters A, Ve, u^o and Dcfr- Since the overall volume 
flow in a closed cell is zero so that the mean particle velocity, (f )p, is equal to fg, the position of 
the center of mass of the heterodyne wing is determined by fg, whereas the osmotic velocity Ueo 
determines its asymmetric broadening due to the solvent flow, and D^q its symmetric broadening 
due to the assumed isotropic diffusion. Finally, A is the frequency-integrated part of the power 
spectrum arising from the positive-frequency super-heterodyne wing. By the sum rule in Eq. (l20l) . 

A = /,(/,(q)), (58) 

so that A is determined by the pref actor of the super-heterodyne wing. 

rV. EXPERIMENTAL RESULTS 

Before analyzing the results of our electrophoretic scattering experiments, we first describe our 
findings for the field- free static structure factor of deionized systems. Next we discuss our elec- 
trophoretic flow experiments starting first with suspensions of non-interacting colloidal spheres 
and progressing then to fluid-like systems of strongly interacting spheres. Finally, crystalline sys- 
tems will be shortly addressed. 

A. Static light scattering 

A standard static light scattering experiment was employed to obtain the measurable static 
structure factor, Sm{(i), from the division of the (background corrected) scattered mean intensity 
(Isil)) —h{<l)u20 of an ordered suspension by the (background corrected) scattered intensity from 
a non-interacting suspension of the same species (at salt concentration of c = 4 x 10"'^ mol/1) 
weighted with the dilution ratio rinon-interactmg /^ordered- The background intensity I{q)B20 of par- 
asitic stray light was determined from measurements of I(q) on the pure solvent. The static 
structure factor shown in Fig. [5] is for a deionized suspension of PnBAPS68. Experimentally, 
the pronounced principal peak and the consecutive secondary oscillation are resolved. From a 
comparison of the signal obtained for the salty suspension to a Rayleigh-Debye-Gans form fac- 
tor we found an angle-independent background intensity of 0.15 in relative units due to multiple 
scattering. This offset was subtracted from the static structure factor measured on the deion- 
ized suspension. A fit of the experimental structure factor using the decoupling approximation of 
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FIG. 5: Experimentally determined static structure factor Sexpil) of a slightly polydisperse, deionized 
suspension, with cr = 68 nm, = 8.38 x 10^^, and s = 0.05, in comparison to the decoupling approxima- 
tion structure factor Soiq), and the structure factor Sid{q) of an ideally monodisperse system calculated in 
RMSA. The effective particle charge has been adjusted to match SD{q) to the peak height of the experimen- 
tal Sexpiq), by adjusting the effective charge to 575 e. Note that X{q) « 9s^ is practically g-independent 
in the experimental g-range. 

Eq. (|48l) for a polydispersity s = 0.05, and an ideally one-component S{q) obtained from the 
rescaled mean-spherical approximation (RMSA) and a DLVO-type screened Coulomb pair poten- 
tial, describes the multiple- scattering-corrected data overall quite well. Note that the value of the 
scattering wave number, q*, in our super-heterodyne flow experiments is in the low-g region of the 
structure factor, with q* /qm ~ 0.4. From the theoretical fit, we have extracted the low-g values of 
Soi^q = 0) ~ 0.03, and S{q = 0) ~ 0.0034, respectively, showing that the incoherent contribution 
dominates S£,{q, r) at very small q. A similar theoretical result is found for the silica particles. In 
particular for PnBAPS68 the experimental data are above the theoretical expectation, indicating 
that there is a larger incoherent contribution. Such a contribution is frequently observed in static 
light scattering and presumably results from optical inhomogeneiety of the particles. 

B. Electrokinetic flow of non-interacting particles 

At a number density of n = 3.95 x 10^^ (0 3.3 x 10""^) and a salt concentration of c = 
5.15 /xmol/l corresponding to a reduced screening parameter contribution Kga ^ 0.97 a suspension 
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of Si251 spheres is essentially non-interacting. The measurable static structure factor, Sm{(i), in 
this system is practically equal to one. Fig. [6^ shows a series of measured power spectra taken for 
increasing field strength. The experimentally accessible frequency range is restricted to values in 
the vicinity of the Bragg frequency of 2 kHz, whose value corresponds to a zero Doppler shift for 
the super-heterodyne part. Here and in the following figures, the frequency / is taken relative to 
the Bragg shift, i.e., f ^ f — A/^, so that only the Doppler- shifted frequency is displayed. Near 
the frequency / = corresponding to the Bragg shifted value, the data in Fig. [6] are more noisy 
and some spectra show additional intensity. This feature is not reproducible and an experimental 
artefact of still unknown origin. Possibly it stems from stray light reflections containing some low 
frequency contaminations superimposing with the reference beam. The overall spectral shape, 
however, clearly reflects the diffusion-convoluted velocity distribution of the heterodyne wing 
depicted in Fig. |4] and based on Eq. (l57l) . There is a linear shift of the center of mass of the 
spectra with increasing field strength E, corresponding to an increasing electrophoretic velocity 
\ve\ = /ie E. The width of the spectra also broadens, as the electro-osmotic velocity, Ueo, near 
the cell wall increases linearly with the field strength. In addition, there also is a symmetric 
broadening of the spectra due to isotropic diffusion. These features have been anticipated also 
in earlier heterodyne dynamic light scattering results [|l8n . but there they have been masked by 
homodyne contributions and overlap with the complementary 'left-handed' heterodyne wing nloA . 
With our super-heterodyne setup, these features can be clearly discriminated, allowing thus a fully 
quantitative evaluation. 

Quite interestingly, we find experimentally that the shape of the spectra of non-interacting 
spheres is given by a single master curve, in accord with the theoretical expectation when Des oc 
E,v (X E and P{v) oc l/E. This can be seen in Fig. [6b where the spectra superimpose on each 
other when the scaled intensity C^^|^^{q, f)xEis plotted versus the scaled frequency f /E. Except 
for the spurious contributions of low frequency noise near the origin, all spectra coincide nicely 
within the experimental errors. This shows that the simple superposition of velocity-weighted 
single-particle Lorentzians used in Eq. (l57l) works quite well for the considered systems of non- 
correlated particles, where 6*^/(5) ^1. In other words the measured spectra can be directly 
interpreted as velocity distributions broadened by diffusion-like mechanisms characterized by D^n- 
This allows us to include a velocity scale, at the top of Fig. [6^ by multiplying / with 2tx / q = 0. 796. 

We have fitted all measured spectra by Eq. (l57l) using the four independently adjustable param- 
eters: Ve, Ueo, and integrated scattering signal strength A. The fits describe the experimental 
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FIG. 6: Power spectra obtained for non-interacting Si251 spiieres at a salt concentration of c = 5.15 ^mol/1 
and a number density of n = 3.95 x 10^^ m^'^. The frequency is taken relative to the Bragg-shift A/b- 
The smooth solid lines are fits of Eq. [57] to the data. The power spectra are well described in all cases, thus 
the superposition procedure works well. The upper scale shows the corresponding velocities obtained from 
V = 2tt f /q. b: scaled spectral power versus reduced frequency. All spectra nicely coincide, showing that 
in all cases the flow is parabolic. 

power spectra of non-interacting particles quite well. This implies, as expected, that the flow pro- 
file in these systems is of parabolic shape. The results for the optimized parameters of these fits 
are included in Figs. |7^-c and compared with the results for suspensions of strongly interacting 
particles that will be discussed in the following subsection. Both the electrophoretic and electro- 
osmotic velocities deduced from the fit increase linearly with increasing field strength (see Fig. 
|7^). This is an important first check, showing that we are in the linear response limit of electroki- 
netics and that the applied field strength is much smaller than the field within a particle double 
layer. Our fits further yield Ueo > \ve\, corresponding to a strong electro-osmotic shear flow close 
to the wall. 

Furthermore, a roughly field-independent effective diffusion coefficient is obtained for the non- 
interacting case with Dcs = 4 x 10^^^ m^/s. The coefficient has a value somewhat larger than 
the Stokes-Einstein value Dq ^ 1.7 x 10^^^ m^/s of the Si251 spheres, and shows the indication 
of a slight increase with the applied field. A more strongly pronounced and linear increase of 
Deff with the applied field has been observed in an earlier study on low salt, but non-interacting 



polystyrene latex suspensions UM, where -Deft was found to extrapolate approximately to -Do in the 
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FIG. 7: Fit parameter results in dependence on the applied field strength, a: electrokinetic velocities, 
|fe| (down triangles) and Ueo (up triangles) for non-interacting Si251 at n = 3.95 x 10^^ m^'^ and c = 
1.32 /imol/1; b: comparison of effective diffusion constants, Df.ff, for the following suspensions: Si251 
at n = 3.95 x 10^^ and c = 5.15 ^umol/l (up triangles), Si251 at n = 3.95 x 10^^ m^^ and 
c = 1.32 yumol/l (down triangles), Si251 at n = 3.95 x 10^^ and c = 1.07 /xmol/1 (diamonds) 
and PnBAPS68 at n = 0.47 fim^^ and c = 0.2 |Umol/l (closed circles). The two horizontal lines indicate 
the Stokes-Einstein values Dq = 6.30 x lO^^^^Vs and Do = 1.71 x W'^^m^/s for PnBAPS68 (dotted) 
and Si251 (dashed), respectively, c: comparison of integrated spectral intensities A for Si251. Symbols as 
before. The principal peak of S(q) in the three suspensions takes values of about 1,1.3 and 1.7, respectively, 
indicating increasing particle correlations. 

zero-field limit. The enlarged value for D^s at small E may have several origins. Most probable, 
transient time broadening, and/or an increased sample temperature due to Joule heating expected 
at this elevated salt concentration cause the enlarged Z^gff- The possible additional increase with 
increasing field will be discussed in detail below. 
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The frequency-integrated power spectrum of the non-interacting system is found to be essen- 
tially field-independent (up triangles in Fig. U^). According to Eq. (l58l) . the frequency-integrated 
intensity of the super-heterodyne wing for non-interacting monodisperse particles is given by 
A (X n P{q), where P(q) is the single-sphere form factor. It is independent of the (normalized) par- 
ticle velocity distribution and Deft- An enlarged E stretches the super-heterodyne wing due to the 
increased Ueo (solvent flow profile becomes more parabolically stretched out), and frequency- shifts 
its center of mass due to the increased value of Ve, but the frequency-integrated super-heterodyne 
part, i.e., the 'area' A under the wing, should remain constant unless P{q) or n would change 
significantly in the presence of a field. 

C. Electrokinetic flow of fluid-like suspensions of interacting particles 

We now proceed from non-interacting suspensions to systems of increasingly pronounced 
short-range order and decreasing isothermal osmotic compressibility. Upon reducing the salt con- 
centration at constant n = 3.95 x 10^^ m~'^ for Si251, or increasing the particle concentration of 
deionized (c = 0.2 yumol/l) PnBAPS68 spheres from n = 0.47 yum~^ to n = 4.7/im~^, the peak 
heights of the field-free static structure factor have values of 1.3 and 1.7 for Si251, and 2.0 and 2.6 
for PnBAPS68. Thus all these samples are in the isotropic fluid state with their short range order 
becoming more and more pronounced in that order. We first discuss the spectral shapes and then 
turn to the effective diffusivities and the frequency-integrated power spectrum. 

Let us first make a general remark on spectra recorded at moderately strong interactions. Ac- 
cording to the sum rule in Eq. dSB, fi'i?(q, 1") ~ |c|^A^/^(g) [9s2G(q, r) + ^(q, r)] ,theIACFis 
a weighted superposition of coherent and an incoherent part. An expression for the coherent part 
in an external field is not yet known. Correspondingly, the integrated spectral power in Eq. (l58l) 
is given by A (x nP{q) [9s^ + ^(q)] for slightly size-polydisperse systems and, more generally, 
hy A (X np{q) Sm{^)- With increasing interactions the low-g form of S'(q) decreases but the the 
incoherent scattering contribution may not yet dominate the coherent one. In such a situation, the 
major assumption underlying Eq. (|52|) . namely the dominance of the incoherent scattering part, 
does not necessarily apply. Rather, we expect also the coherent (steady-state) dynamic structure 
factor, ^(q, r), to contribute to the spectrum. And since to date no exact theoretical expression is 
known for this quantity and its Fourier transform, one should be cautious in interpreting spectra 
of systems of moderately interacting particles as velocity distributions. On the other hand, any 
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FIG. 8: Suspensions of weakly and moderately strongly correlated particles of Si251 at n = 3.95 x 
-|^q16 jj^-s ^. field-strength scaled spectra at c = 1.32 ;umol/l, for field strengths as indicated; b: as in 
a but for c = 1.07 /imol/1; c: comparison of the E = 68.3 V/cm spectrum in b: with the fit curve derived 
from Eq. ([57]). Note the extra spectral power in the experimental data at moderately small frequencies. 

(field-dependent) deviation from the expected spectral shape and amplitude may indicate interest- 
ing underlying physics. 

In Fig. [8^, we show the field-scaled power spectra for the weakly interacting silica spheres. 
Again all spectra superimpose neatly. Only minor changes of the spectral shape are observed, 
mainly appearing as a slight increase in the peak region as the field is increased. In Fig. [8h, 
five spectra are superimposed for the moderately strongly interacting silica suspension. The field- 
dependent changes of the spectral shape in Fig. 8b are more pronounced than before. In particular 
one observes a relative increase of spectral power for the high-shear, low velocity part of the spec- 
tra. While this feature is essentially absent at the lowest fields and high-quality fits are obtained, 
the fits at large fields are considerably worse. This is highlighted in Fig. [8]:, where we compare the 
(unsealed) experimental super-heterodyne signal to a fit of Eq. (l57l) . The difference between ex- 
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pected and measured spectral power distribution in the interval —200 Hz < / < 200 Hz is clearly 
visible. Neglecting this region, a reasonably accurate description of the data is still achieved, com- 
promising between peak intensity and diffusive-like broadening. A fit by Eq. (l57l) used for is 
based on the assumption of homogeneously distributed, identical scatterers and a parabolic flow 
profile. Thus either local changes of scattered intensity or deviations from parabolic flow may 
cause the deviations observed. The restricted frequency range for the relative enhancement of the 
experimental spectral power points to a local effect. One possible cause could be shear- or wall- 
induced changes in S{q, r). In principle, it should be even possible to discriminate between local 
effects on the scattering power and global changes of the flow profile. To this end we checked 
for the linearity of the increase in the velocity corresponding to the center of mass of the spectra 
with increasing field strength. An increased scattering power at low frequencies would indicate 
a decreased apparent velocity. In contrast to this, a changed flow profile would no be such an 
indication, since for closed cell boundary conditions the integrated solvent motion must average 
to zero. In the present case, the observed deviations from a linear increase of the velocity with 
increasing field was still below the experimental uncertainty in Ve. 

In Figs. |9^-b, we show the scaled spectra for low-salinity suspensions of strongly interacting 
PnBAPS68 spheres. The extension of the spectra in terms of scaled frequency is more narrow 
than in the spectra considered before, and they start at about zero frequency. The electro-osmotic 
velocity in this case is much smaller, and on the order of the particle electrophoretic velocity. In 
turn, this implies a less curved parabolic solvent profile stretching to smaller velocities in the cell 
center. Interestingly, the low-/ overshoot of the experimental power spectrum is absent in these 
samples. Second, in particular for the more concentrated sample, the overall spectral power in- 
creases significantly with the applied field. Still, all these spectra can be fitted rather well using Eq. 
(|57]) . and Fig. 9c gives an example of the excellent fit quality. The simple superposition procedure 
in Eq. (|57l) again works well and the parapolic flow profile is confirmed. Here the reduction of 
the low-g form of 5'(q) is sufficient to obtain a super-heterodyne spectrum dominated by inco- 
herently scattered light, in accord with the requirement underlying Eq. (|52l) . Thus, like for the 
non-interacting case, we again may interpret these spectra as velocity distributions. This main re- 
sult of the experimental part of our paper on flow measurements is corroborated by the observation 
that the experimentally determined electrophoretic velocities are in excellent agreement with both 
state-of-the-art analytical calculations [|57n . and combined Monte Carlo and Lattice-Boltzmann 



computer simulations on the Primitive Model level [|53, 
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FIG. 9: Fluid-like suspensions of strongly interacting PnBAPS68 particles, a: scaled spectra for PnBAPS68 
at n = 0.47 ^um^^ and c = 0.2 ^mol/1, for field strengths as indicated; b: scaled spectra at n = 4.7 /um~^ 
and c = 0.2/imol/l for field strengths as indicated, c: background corrected spectrum for n = 4.7 ^um"^ 
and E = 35.1V /cm, with fit curve according to Eq. (l57l) d: field-strength dependence of the integrated 
spectral power for different particle concentrations as indicated. 

Given this, a further comparison of the other two fit parameters reveals some interesting de- 
viations which seem to contain information on interesting underlying physics, not yet fully un- 
derstood. The effective diffusion coefficients. Den, of the silica systems at different salinities and 
as a function of increased field were compared in Fig. |7}5 to those of PnBAPS at n = 0.47m~^. 
For the three Si251 samples, the zero-field extrapolated values are above the Stokes-Einstein re- 
sult (dashed line) with a mild decrease with increasing particle interactions. It is well known that 
the (equilibrium) long-time self-diffusion coefficient decreases below the Stokes-Einstein value 
with increasing interactions, whereas the opposite trend appears for the collective diffusion co- 
efficient Dc = Dq H(0)/S{0) yj]. Our finding therefore nicely illustrates the sum rule in Eq. 
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FIG. 10: Integrated spectral power for a shear melting suspension of PnBAPS68 at n = 5.8 ^m^'^. A field 
of = 46.1 V/cm was applied at t = 25 s. After the initial jump from its zero-field value, A stays constant 
irrespective of the system structure 

(|5TI) . as it shows an originally dominating, but decreasing influence of the coherently scattered 
light with decreasing isothermal compressibility, which in turn lowers the contribution of collec- 
tive diffusion. For the deionized suspension of PnBAPS68 spheres, linear extrapolation gives a 
value Df.ff{E = 0) ~ 1 x 10~^^m^/s well below Dq (dotted line), and in accord with the expecta- 
tion for the value of the long-time self-diffusion coefficient. Thus, estimates of the latter quantity 
may be obtained from extrapolation of Deff(E = 0) to zero field, but only for strongly interacting 
fluid systems, where incoherent scattering dominates and additional broadening mechanisms are 
absent. 

Fig. 7b further reveals an increase of the effective diffusion coefficient, D^-ff, with increasing 
field strength. The effect is only slightly visible at low n and large salt concentrations, but be- 
comes more pronounced for decreased salinity in Si251, and increased particle number density in 
PnBAPS68. The field dependence thus increases with increasing interaction strength. Conversely, 
the approximate theory embodied in Eq. (l57l) predicts a field strength independent diffusivity, 
which in fact was observed in the non-interacting case. Li the deionized case. Joule heating cannot 
be responsible for this dependence, as it vanishes with decreasing salt concentration. Transient 
time broadening, which couples the particle velocity to a non-diffusive peak-broadening, may in 
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principle play a role, but the experimentally observed broadening is too large for this fl8^. In ad- 
dition, the effect is also pronouncedly present for PnBAPS68 at strong fluid-like order, where the 
electro-osmotic profile, and thus the absolute velocities are much smaller. We therefore exclude 
both mechanisms from causing the observed field dependence of the effective diffusion coefficient. 
A number of possible causes will be discussed at the end of this sub-section. 

The field dependence of the integrated spectral power is shown in Fig. 7c for the three silica 
systems and in Fig. 9d for four differently concentrated suspensions of PnBAPS68. The integrated 
power for the silica systems does show roughly field independent values, which differ for each 
suspension. Also for PnBAPS68 each system shows a different overall intensity but here for 
the more concentrated systems a trend of A to increase with the applied field is obvious. The 
first observation is traced back to a detail of the experimental procedure. Each system has a 
different scattering power at low q. This is either due to the different isothermal compressibility 
(for Si251 at different salinity) or different particle concentration (PnBAPS68). With changing 
scattering power of the suspension, also the reference beam intensity has been adjusted to obtain 
an optimized signal-to-noise ratio. This manual adjustment, however, is not precisely defined, 
leading to differences in A for different suspensions. Hence the absolute values for A should not 
be compared. For a given suspension, however, the optical adjustment is kept and hence the field 
dependence can be studied without ambiguity. Therefore the observation of an increased spectral 
power with increased field strength is an interesting finding, which is possibly connected to the 
increase of the effective diffusion coefficient with the applied field. Recalling that in general, 
A (X n p{q) S'a/ (q), a field induced density fluctuation on a suitable length scale will lead to an 
enlarged value of A with increased field. To pursue this idea, further, we performed a measurement 
of the frequency-integrated power spectrum as a function of time for a shear-melting suspension 
of PnBAPS68 dXn = 5.8 jjimT^ with the spectra shown in Fig. 10. The integrated spectral power 
jumped discontinuously from its zero-field value to a 30% larger value upon application of the 
field. The fully crystalline state was stable for about half a minute, then the crystal slowly melted 
inwards and the final stationary state of a shear-molten structure in the mid-plane of the cell was 
reached after some five minutes. During the whole process, A stayed practically constant. This 
suggests that the increase of spectral power observed in the field-dependent measurements is not 
due to a change in the structure. Rather, the initial jump seems to indicate a correlation with the 
presence of the electric field. 

Several previous studies were concerned with the coupling of diffusivity and/or density flue- 
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tuations to external fields. Qiu et al. n64\ studied the long-time self-diffusion of non- or weakly 
interacting charged spheres under low frequency, oscillatory Couette flow. With increasing shear, 
they observed anisotropic self-diffusion, with the mean-squared displacement in the direction of 
the flow field growing faster than linear in time. This was attributed to Taylor dispersion driven 
by the position-dependent flow [6^. In addition, a non-convective enhancement of diffusion was 
found also perpendicular to the flow direction, characterized by a diffusion constant that increases 
with the shear rate. This shear-induced enhancement of self-diffusion in strongly sheared, di- 
lute suspensions of charged particles was attributed by Breedveld and Levine to the electrostatic 
breaking of the time reversal symmetry of particle trajectories in low Reynolds number flow [|67|l . 
Stokesian dynamics simulations for non-dilute, sheared suspensions of hard spheres also show a 
pronounced shear-induced enhancement of long-time self-diffusion perpendicular to the flow di- 
rection 116 811 . Conversely, recent analytic calculations and Brownian dynamics simulations (without 
hydrodynamic interactions) of collective diffusion of hard spheres ll42ll found the shear-induced 
change in the collective diffusivity tensor away from isotropy to be quite small at small shear rates 
(Peclet numbers). Not much is as yet known about suspensions of strongly interacting particles 
under shear. But, due to the overall parabolic flow profile for fluid-ordered suspensions in the 
present experiments, we expect an inhomogeneous spatial distribution of observable effects. 

In a recent study on the electrokinetics of charged spheres in a quiescent solvent, Araki and 
Tanaka ll52ll performed numerical calculations based on a so-called fluid particle dynamics (FPD) 
method, where the colloidal particles are described as highly viscous liquid spheres and the sol- 
vent and microions are treated as continuous fields akin to standard electrokinetic theory. The FPD 
method accounts for many-particle hydrodynamic interactions. Within the linear response regime 
of the field-dependence of the electrokinetic velocities, the authors found hydrodynamically in- 
duced velocity fluctuations which increased with the field strength. Clearly, such fluctuations will 
directly mimic an enhanced diffusivity. For low field strengths, the velocity fluctuations decreased 
monotonically with increasing salinity but at larger field strength the scaled velocity fluctuations 
pass through a minimum at some intermediate salinity value. The authors further pointed out that 
these fluctuations significantly alter the pair correlation function, i.e. they translated into density 
fluctuations. It will be interesting to explore how such fluctuations, which should likewise grow 
with the electric field, behave on large spatial scales and translate into alterations of the low-g 
S{q). 

Finally, an enlarged apparent diffusivity and corresponding global density fluctuations of non- 
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hydrodynamic origin, can be expected when fluctuations in relative particle velocities are caused 
by a non-negligible mobility polydispersity. Here particles move at different constant speeds which 
introduces a broadening of the velocity distribution. Similar to hydrodynamic broadening, this 
effect is independent of the absolute particle velocity but proportional to the relative velocity be- 
tween neighboring particles i and j which will increase with the applied field: At> = (/Xj — fXj)E. 
In the former case the effective diffusivity is enhanced, in the latter the velocity distribution is con- 
voluted with the mobility distribution (e.g., of Gaussian shape for a narrow distribution). As both 
types of independent spectral broadening mechanisms cannot be discriminated properly within 
the spectra, polydispersity of mobilities may mimic an increased diffusion. On the other hand, 
also the suspension microstructure is likely to be affected by the systematic differential motion of 
neighboring particles. If density fluctuations occur at the appropriate scale, the relative amplitude 
of the coherent scattering contribution would be enhanced. This spatially global effect should 
become stronger with increasing field strength, but could also decrease with increasing particle 
interactions, as the isothermal osmotic compressibility decreases and the structural ordering be- 
comes more pronounced. A discrimination between all these possibilities may be attempted by 
applying our method to different flow situations, with particular focus on those, where the external 
field couples differently to particles, micro-ions and solvent. 



D. Electrokinetic flow of crystalline systems 

In this last subsection, we shortly address crystalline systems as formed in deionized suspen- 
sions of PnBAPS68. In this case, the finite shear moduli and the different viscosities of the dif- 
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ferent phases coexisting under shear cause deviations from a parabolic flow profile [131 , 
We therefore expect that Eq. (|57l) cannot be used. In fact, the shape of the power spectrum in 
the crystalline state at low electric field strengths is quite different to that of the stationary state 
at large fields, where the systems shear melts within a few minutes. This is seen for the lowest 
applied field strength in Fig. [TT^ . for PnBAPS68 at a particle density ofn = 5.8 jim^^. A qual- 
itative estimate of the flow behavior was obtained from microscopy. At low field strength E, one 
observes some kind of double plug flow, lubricated by a layer crystal region. At a value E = 22.05 
V/cm, the layer crystal fills the complete cell and the flow profile takes a parabola-like shape. At 
E = 34.6 V/cm and larger, the system is shear-molten at mid-cell height. The system of Fig. [TOb 
is at n = 26.5 yum~^ and stays fully crystalline even at the largest applied field strength. Micro- 



40 



0.210 r 



PnBAPS68 n = 5.8^lm 
57.60 V/cm 




4 8 12 16 
(f / E) / Hz cm V ' 



20 



1.6 

o 

> 1.2 

CO 

lu 0.8 

X 

^1 0.4 
o" 

0.0 



b 


mil iiiiirn-^ 


PnBAPSeS 
n = 26.5^Jm ' 

\ 11 9.2 V/cm 

'■WI...II 


- 


-J 


^ 83.4 V/cm 




59.6 V/cm 


_j 1 


35.9 V/cm 


1^ J\ 11.9 V/cm 



4 8 12 
(f/E) / HzcmV ' 



16 



20 



FIG. 11: Scaled power spectra for two equilibrium crystalline suspensions of PnBAPS68. Spectra are 
shifted upward for clarity, a: n = 5.8 ^m~^. With increasing field strength this system consecutively shows 
double plug flow, layer crystal plastic flow and complete shear melting; h: n = 5.8 ij,m^^. With increasing 
field strength this system shows double plug flow but no shear melting. All flow profiles, even those in the 
shear molten state, deviate strongly from the parabolic flow profile entering into Eq. (|57] ) . 

scopic observation here reveals a double plug flow in opposing directions. The directly observed 
flow behavior is mirrored in the spectral shapes, and none of these spectra can be fitted well by a 
parabolic flow profile according to Eq. (l57l) . Interestingly, this applied also to the systems which 
were completely shear-molten at mid-cell height. Height dependent measurements show that this 
is due to an incomplete shear melting in y-direction, altering thus the aspect. 

Clearly, the homogeneity assumption used in the derivation of the scattering theory is valid only 
locally in (partially) crystalline systems, and invalidated at each slip line and/or phase boundary, 
while at the same time the long time self-diffusion coefficients are expected to differ considerably 



for different suspension structures ll69ll . Therefore, Eq. (|57l) has to be applied with great caution in 
both integral and local measurements. For the low diffusivity case of crystalline or partially shear- 
molten states under hydrostatically induced flow through a capillary, the comparison of local and 
integral spectra, however, suggests that in a qualitative way the latter can still be interpreted as 



velocity distributions and estimates, e.g., of the electrokinetic velocities be obtained [j33l]. More- 
over, Fig. 10 has shown that the spectral power does not change upon changing the suspension 
structure. This allows further quantitative evaluation of flow profiles, as the intensity at a given 
frequency is then proportional to the number of scatterers moving with the corresponding veloc- 
ity. Hence, within the uncertainty of unknown but presumably small diffusional broadening, quite 
complex flow profiles should be analyzable at least on a qualitative level. In fact, flow profiles con- 
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sistent with theoretical expectations for multiphase flow have been extracted from point- by-p oint 



measurements of local velocities in pipe flows under hydrostatic pressure differences |128, 



V. CONCLUSIONS 

Integral heterodyne velocimetry is an established method to measure flow in colloidal systems. 
It was applied in the past both in electrokinetic and rheologic contexts to study electrophoretic 
mobilities and multi-phase flow, and the theoretical expressions derived for non-interacting sus- 
pensions were applied to estimate velocity distributions and flow profiles. In the present paper, 
we went beyond this ad-hoc ansatz and extended the experimental method to super-heterodyning. 
Moreover, we could complement this by a derivation of a theoretical expression accounting for 
structure formation and the unavoidable polydispersity of real colloidal systems and a careful dis- 
cussion of the performance and scope of this approximate approach. Our analysis shows that for 
small size polydispersity and sufficiently low isothermal compressibility (more precisely: low-g 
coherent scattering intensity) of the suspension, the signal in the low-g limit should be dominated 
by incoherent scattering. Experimentally, a low compressibility was realized by using strongly 
interacting low-salt suspensions, in the fluid or crystalline state. In our theoretical considerations, 
we first assumed a constant particle electrophoretic velocity and a homogeneous suspension. Ex- 
perimentally, we studied electrokinetic flow in closed cells yielding parabolic flow profiles. Spatial 
variations of the velocity were theoretically accounted for by convoluting a theoretical expression 
for the particle velocity distribution with the single- velocity spectral power form. In this way, we 
obtained an excellent description of the measured power spectra for fluid-ordered systems. 

Under the assumptions made in the scattering theory presented here, the derived expression for 
the incoherent spectra is identical in form, but not in value, to that of the non-interacting case. It 
can be fitted to experimental spectra using four independent parameters of which two (effective 
diffusion coefficient and integrated spectral power) are disregarded in the further evaluation for the 
electrokinetic velocities. This neglect works well even in the case of crystals and of field-induced 
structural changes [26, 27, 28]. The electrophoretic mobility for shear-melting systems was found 
to be constant, irrespective of the applied field strength and the stationary system microstructure. 
This finding is fully consistent with recent simulations including hydrodynamic interactions [|54|| . 
So there already is a practical assessment of this approach for the complicated case of electroki- 
netic flow and shear. Inclusion of different experimental situations, like hydrostatic tube flow. 
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Couette-shear, sedimentation, and even shear-banding flows appears to be rather straightforward. 
Hence the presented method promises to become a valuable tool in the studies of colloidal disper- 
sions under flow. 

A deeper theoretical understanding, however, remains as a challenge, to be addressed in future 
work. In particular, the observed dependencies of the effective diffusion coefficient and the spec- 
tral power underline the importance of fluctuations, which have not been considered so far in the 
above mentioned evaluations but should be included in our simple theoretical model. Irrespective 
of their origin (shear-induced Taylor dispersion, hydrodynamic coupling by micro-ion currents, 
polydispersity in the electrokinetic mobility or sedimentation velocity), it is reasonable to expect 
that differential motion of neighboring particles will affect the suspension micro structure and may 
couple to density fluctuations on the scale of nearest neighbors distances and larger. So in all these 
cases, we expect an enhanced contribution of coherent scattering to the spectra. In particular, an 
extension of the theoretical treatment to the coherent dynamic structure factor, 5'(q, r), is highly 
desired, which also includes a worked-out calculation of its self-part ^(q, r). This should also 
include strongly interacting suspensions with fully developed fluid microstructure and polydisper- 
sity effects. This, of course is quite demanding, since there is no general expression available to 
date for the steady-state dynamic structure factor S'(q, t) in driven systems covering all cases of 
driving fields ranging from simple shear over sedimentation, where we have a relative motion of 
colloidal particles against the backflowing solvent, to electrokinetic motion which combines shear 
flow with particle motion against the solvent. In all three cases the coupling between structure 
and diffusivity via direct interactions and hydrodynamic interactions is different, and in any case 
very complex. A practical aid in this task could again be given by the experiment, by comparing 
angle-dependent measurements in the three different situations. This could clarify whether the 
changes in the power spectrum are caused by the particle motion relative to the solvent, the shear 
flow or the applied electric field. Special interest may further be given to the flow forms in binary 
mixtures of low-polydispersity samples. By a suitable choice of particle material, size and charge, 
each of the mentioned external perturbations may act differently on the two species and induce 
clearly separated signatures in the spectra. 

The present paper has reported on an important methodical development, both from the ex- 
perimental and theoretical points of view. It sets the stage for progressing from the quantitative 
characterization of flow in terms of transport coefficients to a deepened and more general un- 
derstanding of the complex interplay between particle interactions, microstructure and dynamics 
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under flow. 
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